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Abstract 

The present study is focused on the superscaling behavior of electron-nucleus cross sections in 
the region lying above the quasielastic peak, especially the region dominated by electroexcitation 
of the A. Non-quasielastic cross sections are obtained from all available high-quality data for 12 C 
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by subtracting effective quasielastic cross sections based on the superscaling hypothesis. These 
residuals are then compared with results obtained within a scaling-based extension of the rela- 
tivistic Fermi gas model, including an investigation of violations of scaling of the first kind in the 
region above the quasielastic peak. A way potentially to isolate effects related to meson-exchange 
currents by subtracting both impulsive quasielastic and impulsive inelastic contributions from the 
experimental cross sections is also presented. 

PACS numbers: 25.30.Fj, 24.10.Jv, 13.60.Hb 
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I. INTRODUCTION 



properties of electron-nucleus scatter- 



In recent years scaling [l[ and superscaling {3], Q] 
ing have been studied in great detail. A first line of investigation has been focused on the 
behavior of experimental data and on the construction from them of suitable phenomeno- 

A second line, developed in 



logical models for lepton-nucleus scattering 
parallel to the first, has instead been focused on more theoretical analyses; namely, the 
superscaling properties of cross sections obtained within specific nuclear models have been 



analyzed with the goals of testing the range of validity of the supersca 
of finding and explaining possible scaling violations j^.Tlol. 
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epton-nucleus scattering in the region of the A resonance has been recently studied in 



2fl 



21] and an extension of the scaling formalism to neutral current neutrino processes has 



also been proposed 
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The general procedure adopted in scaling analyses consists of dividing the experimental 
cross sections or separated response functions by an appropriate single-nucleon cross section, 
containing contributions from protons and neutrons, in order to obtain a reduced cross 
section which is then plotted as a function of an appropriate variable, itself a function of the 
energy and momentum transfer. If the result does not depend on the momentum transfer, 
we say that scaling of the 1st kind occurs. If, additionally, the reduced cross section has 
no dependence on the nuclear species, one has scaling of the 2nd kind. The simultaneous 
occurrence of scaling of both kinds is called superscaling. 

The superscaling properties of electron-nucleus scattering data in the quasielastic (QE) 
region have been extensively studied in jjj, 4] and in [5]: scaling of the 1st kind was found 
to be reasonably well respected at excitation energies below the QE peak, whereas scaling 
of 2nd kind is excellent in the same region. At energies above the QE peak both scaling 
of the 1st and, to a lesser extent, 2nd kinds were shown to be violated because of the 
important contributions introduced by effects beyond the impulse approximation: inelastic 
scattering j^, [(], 25^ . correlations and meson-exchange currents (MEC) in both the lp-lh 



and 2p-2h sectors 
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301 ]. which mostly reside in the transverse channel. 



The variety and complexity of contributions that are present above the QE peak make it 
difficult to analyze inelastic data directly in terms of inelastic scaling variables and functions. 
Any analysis of this type requires some kind of theoretical assumption which allows one to 
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focus on a specific kinematic region, having removed contributions from other processes (to 
the degree that one can). In 7j|, the scaling analysis of electron scattering data was extended 
to the A resonance region. A non-quasielastic 1 (non-QE) cross section for the excitation 
region in which the A plays a major role was obtained by subtracting QE-equivalent (see 
below) cross sections from the data and was found to scale reasonably well up to the peak. 
Phenomenologically determined QE and non-QE scaling functions were then used to obtain 
predictions for neutrino cross sections at similar kinematics {?], Q] • This approach has been 
referred to as the SuperScaling Analysis (SuSA). 

In this paper one of our goals is to investigate superscaling, and its violations, in the region 
above the QE peak, starting from the idea presented in |7|. To this purpose, in Sec. HT1 we 
begin by reviewing the basic formalism for scaling studies in the QE region; specifically, 
we summarize the essential features of the so-called SSM-QE model (to be defined in that 
section). We continue in that section by also considering the A region, reviewing and 
extending the SuSA approach of [3]. All available high-quality data for 12 C are reconsidered 
and analyzed by applying a variety of kinematical cuts to illuminate the origins of the 
scaling violations that are observed. We then proceed to a deeper investigation of these 
scaling-violating contributions in the region between the QE and A peaks. In order to do 
so, in Sec. IHII we present a model for inelastic electron-nucleus scattering within the impulse 
approximation based on the same superscaling ideas of p|, extending an earlier superscaling- 
based model for inelastic scattering — this is the so-called SSM-inel approach and has 
a variant denoted SSM-A — see that section for specific definitions. These models are 
used in Sec. IIIIBI to compute non-QE superscaling functions and to compare these with the 
experimental data and with the SuSA fit for several choices of kinematics. By subtracting 
theoretical inelastic cross sections from the experimental data, in Sec. IIVI we then use this 
model to isolate the non-impulsive components of the cross section and analyze their behavior 
in terms of 2p-2h MEC contributions obtained in previous studies. Finally, in Sec. [V] we 
summarize our study and draw our conclusions, including some remarks of relevance for 
studies of neutrino reactions with nuclei. 



1 In [7| this residual was called the "Delta" contribution, assuming the A to be dominant. To avoid 
confusion with later discussions where A-dominance is assumed, in the present work we denote the entire 
residual after the quasielastic contribution is removed by "non-QE" . 
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II. FORMALISM AND PREVIOUS RESULTS 



A. Scaling in the QE region: the SSM-QE approach 



Here we present a summary of the relevant formalism for scaling studies in the QE 
region, focusing on the formulae and results which will be used in the rest of our study. We 
denote this the SuperScaling Model for the QE response functions (SSM-QE). Our purpose 
is to illustrate how scaling ideas can be used to motivate the construction of superscaling- 
based models for electron-nucleon cross sections, in the spirit of 6|, |7| where more extensive 
discussions can be found. 

Within the Relativistic Fermi Gas (RFG) model the only parameter characterizing the 
nuclear dynamics is the Fermi momentum hp. In the following we will retain only the 
lowest orders in an expansion in the parameter rjp = kpjm^^ being the mass of the nu- 
cleon. Within this approximation the RFG longitudinal (L) and transverse (T) quasielastic 
response functions, at momentum transfer q and energy transfer u, can be written as 

<t(m)4^w g s, (i) 

where the scaling function is given by 



fRFG{$) = f^FcW = fRFG 

and the scaling variable ip is 



;i-v 2 )0(i-v> 2 



i 



A 



T 



(1 + X)t + kJt(1 + t 



(2) 



(3) 



with £f 



1 + rjp — 1. In the formulae above we have introduced the usual dimensionless 



variables: k = q/2m N , A = u/2m N and r = k 2 — A 2 . Retaining terms only up to order rjp, 
the functions G^t are given by [5( 



qQE 



JL{z[(i + r)wg p E - 

G% E = - [zW?f + AWgf } 



rr i,p 



+ N[(l + T)Wgf-Wf*\} 



(4) 
(5) 



W(i 2 ) pin) being the single-proton (-neutron) electromagnetic structure functions, which are 



(l,2),p(n) 

given in terms of electromagnetic form factors by 



W l,p(n) 
VV 2,p(n) 



T G 2 M,p{n) ( T ) 



^E,p(n)( r ) + T ^A/,p(n)( r ) 



(6) 
(7) 



v L RT + v T Rr , (8) 



Given the response functions, the QE cross section is then obtained as 

da 

de'dQ ~ ° 

where e' is the outgoing electron energy and Q = (9, </>) is the solid angle for the scattering. 
Here om is the Mott cross section and vl,t are the usual kinematic factors. 

The expressions above suggested, instead of using the RFG scaling function in Eq. (j2J) 
that one may work backwards to obtain an experimental scaling function by dividing the 
QE cross sections by the quantity 

S° E = o M \v L GT + v T G Q T E ] (9) 
and then, for use in discussions of 2nd-kind scaling, multiplying the result by k F : 

Separate L and T scaling functions can similarly be obtained as 



R exp 

fj£(ll>,K)=k F -$. (11) 

In our previous analyses of the world (e, e') data we have found that, for large enough 
momentum transfer (q > 2kp), lst-kind scaling works rather well for values of energy transfer 
u below the QE peak value, ujqe- For large values of uj deviations are observed, coming 
from contributions beyond QE scattering, such as inelastic scattering and MEC effects. A 
separate analysis of the longitudinal and transverse channels shows that these deviations 
mainly occur in the transverse response, while the experimental longitudinal reduced cross 
sections scale much better and up to larger values of uj. This suggests that we can use the 
longitudinal QE experimental scaling function obtained in ^,4] to define a phenomenological 
scaling function. 

In particular, assuming that (i) indeed there is a universal superscaling function and that 
(ii) it can be identified with the phenomenological function extracted from the analysis of the 
QE longitudinal response, we can now work backwards and use this superscaling hypothesis 
to predict cross sections. To be more specific, we define the superscaling model for the QE 
response functions (i.e., what we are calling the SSM-QE approach in this work). This 
consists in using Eq. ([Q), but with 

f SSM-QE ^ _ f QE W _ (12) 
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FIG. 1: Phenomenological fits for the superscaling functions f SSM ~Q E (solid line) and fs^s~^ E 
(dot-dashed) versus the appropriate scaling variable. The RFG superscaling function is also shown 
for comparison (dotted line). 

An important step has been taken here: only the longitudinal cross sections are employed 
in defining the phenomenological scaling function. This choice is based on the fact that 
the transverse cross sections can have significant non-QE or non-impulsive contributions, 
for instance, the former from inelastic excitations of the nucleon (importantly the A) and 
the latter from 2p-2h MEC — see the discussions to follow in the present work. However, 
in lowest order these are not very important in the longitudinal cross section, and thus it 
provides the only opportunity to isolate the impulsive contributions to the nuclear response. 

The phenomenological function f SSM ~Q E employed in the present approach is shown in 
Fig. [1] (solid line), where it is compared with the RFG scaling function of Eq. (T5]) (dotted 
line). Also shown is the phenomenological non-QE scaling function, fs'^s~A^ E to be defined 
below in the following subsection (dot-dashed line). Focusing on the phenomenological QE 
scaling function, one sees that it is significantly different from the RFG result: it is about 17% 
lower at the peak and is asymmetric, having a tail that extends to higher u (in the positive 
ip' direction). In fact, subsequent to obtaining the phenomenological results shown in the 
figure 3|, |4( , relativistic mean field theory (RMF) was employed to obtain theoretical scaling 
functions. This approach is especially relevant at high energies where relativistic effects 
are known to be important. These RMF studies yielded essentially the same longitudinal 
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scaling function as the phenomenological model [llj], and the required asymmetric shape 
of the scaling function was obtained theoretically. We shall return below to comment on 
the RMF transverse scaling function. Still later a so-called semi-relativistic approach was 

n 

pursued [12[ , again yielding essentially the same results. More recently, a deceptively simple 
"BCS-inspired" model was developed |3l|, with the same outcome: a peak height that is 
significantly below the RFG result and an asymmetric shape. Within the flexibility in each 
model and the experimental uncertainties one can say that a single longitudinal QE scaling 
function has clearly emerged. 

In passing, we note that, as is usually done in studies of electron scattering in order to 
reproduce the correct position of the QE peak, in the present study we have introduced 
a small energy shift E s hif t . Within the framework of the superscaling formalism outlined 
above, this amounts to considering a "shifted" scaling variable ip', calculated according to 
Eq. ©, but with A -»• A' = A - E shift /2m N and r -> r' = k 2 - A' 2 . The values k F = 228 
MeV/c and E s hift = 20 MeV have been used in all of the calculations for 12 C presented here 
and in the following sections. 

Having found that the longitudinal QE scaling function is universal, whether treated 
phenomenologically or via models for lp-lh knockout reactions, we now discuss the trans- 
verse QE response. In most approaches one finds that once the single-nucleon cross section 
is removed in defining scaling functions as above the longitudinal and transverse answers 
are basically the same, i.e., one has what has been called scaling of the 0th kind with 
/t^O = /lC0O- However, in what is likely the best model employed so far, the RMF ap- 
proach cited above, one finds that Oth-kind scaling is mildly broken for momentum transfers 
in the 1 GeV region with f T (ij/) > For instance, at q = 500 MeV/c (1000 MeV/c) 

the transverse RMF scaling function is 13% (20%) larger at its peak than is the l ong itudinal 
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one. On the other hand, from analyses of lp-lh MEC contributions 
one sees the opposite behavior, namely, the 1-body (impulse approximation) and 2-body 
MEC contributions to the lp-lh response, which must occur coherently and hence can in- 
terfere, in fact do so destructively and therefore a somewhat lower result is found for the 
total transverse scaling function. Neither of these effects is seen in the longitudinal response 
in leading order. Unfortunately, no single model exists where one has adequate relativistic 
content (as in the case of the RMF approach) and has a consistent way to obtain the MEC 
contributions; indeed, the MEC studies cited above could not be attempted on the same 
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footing as the 1-body RMF computations and could only be undertaken using much simpler 
dynamics. 

Accordingly, we have no better option at present than to adopt some working procedure. 
Henceforth we shall assume that Oth-kind scaling is obeyed and thus take /t^O = /lCVO f° r 
the quasielastic response. One should remember, however, that this may not be completely 
true and that the QE transverse response could be either a bit larger or a bit smaller than 
the one obtained under this assumption. In Sec. IIVI where the scaling-based cross sections 
are compared with data, we shall return to discuss these issues in somewhat more detail. 



B. The SuSA approach to scaling in the A region 



We begin by summarizing the essentials of the SuSA approach taken in [/J, where non- 
QE cross sections were obtained from experimental inclusive inelastic electron-nucleus cross 
sections by subtracting QE cross sections given by the SSM-QE procedure described above. 
Namely the following cross sections 

/ da \ non - QE = f da \ exp _ ( da \ SSM ~ QE 

[de'dn) '{de'dQJ [de'dQJ ^ 

were obtained as a first step. In the earlier work it was assumed that A-dominance could be 

invoked. Namely, in analogy with the QE results of previous section, a model in which only 

impulsive contributions proceeding via excitation of an on-shell A was employed. In that 

model the leading-order RFG expressions for the electromagnetic response function can be 

written as 7L a2] : 



Rl T (K,\) = ^-f A (4, A )Gl T 



with / a Oa) = Irfg^a) and 



A - Tp A 



with 



and with 



Pa = 1 + 



111 -At 



At 



Pa 



m A 

m N 



G 



r<A 



K 

At~ 

1 

2k 



-A 



1 + T P \ + l) 



Aw*. 



(14) 
(15) 

(16) 

(17) 
(18) 



In Eqs. ( TT7ffT8l) the single-hadron iV — > A structure functions are 2 
< = \ (M + I) 2 (2rp A + 1 - /i A ) (G 2 M)P + 3G| in ) 

,„A _ /.. .2 (2rp A + l-/iA) / r2 ^2 , 4 T r 2 

w 2 - + : - — \G MjP + 6G En + 4-2-G CiA 

1 -r "v Pa \ Ma / 

where the magnetic, electric and Coulomb form factors are taken to be 

G M ,p = 2.97# a (t) 
G E , n = -0.03^ A (r) 
G CA = -0.15G A/ , P (r) , 



(19) 
(20) 



(21) 
(22) 
(23) 

- 7 =L= j ■ (24) 

VTT^ (1 + 4.97r) 2 

Starting from these expressions and assuming that the only non-QE contributions arise 
from this A-dominance model one can define a superscaling function in the region of the A 
peak as follows 



with 



9a{t) 



f non ~ QE {%l) L 



da 
de'dVt 



non—QE 



with 



S A 



v L G$ + v T G$ 



We have performed an analysis similar to that presented in [7j and, focusing on scaling of 



17] 



(25) 



(26) 



the 1st kind, we have considered all aval 
cross sections on 12 C 
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able high-quality data of inelastic electron scattering 
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421 ] . The functions f non -Q E we obtain 
are shown in Fig. [2j Note that, as above, we have introduced a small energy shift E s hift- In 
employing Eqs. (TT5l) and ( Ti~6l) we do as in the QE case and replace A by A', and r by r 1 . As 
before, for 12 C the values kp = 228 MeV/c and E s hif t = 20 MeV have been used in all of 
the calculations presented here and below. 

Overall we see a tendency for coalescence below and up to the A peak for some, but 
not all, of the data. Specifically, for kinematics lying below the A peak (ip' A = 0) these 



2 Equations (|19l20p should be taken with A = Z and the p — > A + structure functions and with A — N and 
the n — > A structure functions, and then summed, but since these processes are purely isovector we use 
A = N + Z with one choice for the structure functions. 



9 



non-QE results scale reasonably well given the assumption of A- dominance, showing scaling 
violations at the level of roughly 0.1 units of scaling function, versus the QE peak value 
of about 0.6, namely scaling violations of approximately 15-20%. As discussed in more 
detail later, since we cannot have any inelasticity over much of this kinematic range (being 
below pion production threshold) one must suspect that effects such as from 2p-2h MEC 
contributions are playing a non-trivial role. Nevertheless, accepting this as a measure of the 
potential uncertainty in following the straightforward SuSA approach, in [3] an empirical 
fit to these results, fs^sA > was obtained and then used to predict neutrino-nucleus cross 
sections in the A region. It should be stressed that the assumption of A-dominance is clearly 
only an approximation; in the following sections we present a more microscopic approach 
in which the superscaling approach discussed above for the QE region is extended to the 
inelastic region (denoted the SSM-inel approach; see Sec. lIIIj) and where non-impulsive 2p-2h 
MEC effects are considered separately (see Sec. HV|) . 

Looking in more detail, let us first consider the two bottom panels of Fig. [21 which show 
all available high-quality data for 12 C and the data with momentum transfers q > 500 
MeV/c. We observe that many (but not all) of the data indeed tend to collapse into a single 
function close to the A peak. The spreading of the data is larger than what was observed in 
similar analyses of QE data, but a tendency to cluster (scale) is seen, at least for a subset 
of the data. In order to discuss scaling, and the breaking of it, one may consider three 
different regions. First, in the positive region the spreading of the data increases and 
the data themselves tend to diverge. This behavior is analogous to what happens for the 
QE case for large values of ip' and it is due to the presence of contributions coming from 
higher resonances. Second, for ip' A < — 1 the data form a relatively uniform background 
showing no specific pattern. This is the range where effects from 2p-2h MEC are expected 
to play a significant role (see below). Finally there is the region —1 < ip' A < 0, where the 
spreading of the data is somehow less evident and where both type of scale-breaking effects 
can contribute. 

In a first attempt to disentangle these effects, in the top right-hand and two middle panels 
of the figure we apply a progression of cuts on the data, specifically taking those with 0.5 
GeV/c < q < q cut , where q cut goes from 1 to 2 GeV/c. As the cut tightens we expect to have 
fewer and fewer contributions from higher inelasticities. For completeness, and for a better 
understanding of the whole figure, in the upper left-hand panel we also report the data for 
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FIG. 2: "Experimental" superscaling function f non -Q E for 12 C, obtained by applying the QE- 
subtraction procedure described in the text to the available experimental data for 12 C. The function 
is plotted versus the A scaling variable i// A . Kinematical cuts on the values of the momentum 
transfer q (in GeV/c) are considered, as indicated in each panel. A phenomenological fit of the 
non-QE superscaling function, fs^~j^ E -, is also shown for comparison by the solid line. 

low momentum transfer (q < 0.5 GeV/c). 

The results shown in the different panels seem to indicate that the presence of contribu- 
tions from higher inelasticities corresponds to values of f non -Q E which lie above the average 
scaling function for —0.5 < ip' A < and below it for — 1 < ip' A < —0.5 (for instance, compare 
the top and middle right-hand panels). In particular we observe data sets which seem to 
cross the average function around ip' A = —0.6. They correspond to JLab cross section data 
taken at an incident energy of 4.045 GeV and scattering angles between 23 and 74 degrees, 
for which there is indeed a strong overlap of the A and higher inelastic contributions. 

The observations above suggest that, if we are interested in obtaining a phenomenological 
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SuSA scaling function for the A region alone, fs < ^~A* E , these highly inelastic data sets should 
be excluded from the fit. Such a fit, similar to that obtained in [7|, is indicated in Fig. [2] by 
the solid line, and in Fig. [I] it is compared with the phenomenological fit for the QE region 
and, for reference, with the RFG scaling function. 

We observe that fs < ^sA^ E differs significantly from f SSM ~Q E _ This is expected, because, 
besides incorporating initial-state dynamics, the phenomenological non-QE scaling func- 
tions certainly contain additional effects, such as those due to the finite width of the A 
resonance, as well as potential 2p-2h MEC contributions. However it is interesting to inves- 
tigate whether or not these differences can be explained only in terms of kinematics and of 
trivial effects, such as the finite width of the A, or whether also differences in the nuclear 
dynamics at the QE and A peaks can contribute to them. In order to address this issue we 
need to introduce some model for the cross sections in the A region, and we will present 
this in the next section. 

Let us conclude this section by introducing the phenomenological SuSA model for the A 
region [?J mentioned in the introduction. Following the approach used in previous section 
for the QE case, we can obtain the response functions for A excitation from Eqs. (114j) . by 
substituting the RFG expression for f A with the phenomenological fit obtained from the 
data, namely 



R S ^ A ~ A in, A) = Trfs°uSA QE ^A)Gi T . (27) 

This model was tested in [7| for electron scattering over a range of kinematics, showing 
agreement with the data at the level of 10% or better. 



III. SSM-BASED MODELS FOR THE INELASTIC REGION 

In this section we develop a model for the response functions in the inelastic region lying 
above the QE peak, basing the approach on the assumption of universality of the superscaling 
function, i.e., using the same SSM approach employed for the QE region. This will allow us 
to address two issues. On the one hand we will explore the origin of the difference between 
the phenomenological scaling functions obtained by fitting the data for the QE {^j SSM -Q E ^ 
and A {fsuSA QE ) regions, as discussed in the last section. We will start by assuming that 
this difference can be accounted for only by kinematics and finite width effects, and we will 
compare the scaling function obtained under this hypothesis with the experimental one. 
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On the other hand, we will investigate further the role played by contributions from higher 
resonances in producing the scaling violations shown by the experimental f non ~Q E in the 
region — 1 < ip' A < 0. 

As the model we present here is based on the impulse approximation, it will not allow us 
to investigate MEC effects in the < — 1 region directly, but it will turn out to be useful 
later in Sec. IIVI in presenting the experimental data in a different and more focused way. 



A. Formalism 

We follow closely the approach of {f| , where a microscopic model based on the RFG and 
on superscaling was used to study highly-inelastic electron-nucleus scattering. The RFG 
expressions for the inelastic nuclear response functions can be written as |6j]: 



R 



inel 



f'2 



Kp J ix\ 



inel 



(28) 



where ipx is obtained from Eqs. (|T5l) and (|T6|) for a generic invariant mass Wx of the final 
state reached by the nucleon, namely by replacing /ia with nx = Wx/rrix- The quantities 
G™?, neglecting terms of order rjp and higher as before, are given by 



G 



inel 
L 



m N — {Z [(l + rp|; 



w 



7,P 



W 



+ N (1 + rp\)w n 2 



riinel 



m N -{Zw p 1 + Nw^} 



K 



(29) 
(30) 



where wi^ are the inelastic single-nucleon structure functions, which depend on two vari- 
ables, the four-momentum transfer Q 2 and the invariant mass Wx or, equivalently the 



single-nucleon Bjorken variable x = \Q 2 \/[W 



in 



N 



Q 2 ] (see also 



Note that the 



inelastic structure functions have dimension of E 1 , at variance with the previous QE and 
A cases: for this reason we indicate them as w. The integration limits in Eq. (|28|) are given 
by 



/ii = 1 + ^ 

/i 2 = l + 2A-e 5 



(31) 



with /i„ = m^jmfq and where es = Es/m^ is the dimensionless version of the nucleon 
separation energy. The first limit is simply the threshold for pion production, while the 



second was derived in 
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Following a procedure analogous to that illustrated for QE scattering, we can now gen- 
eralize the RFG by making the substitution 



jSSM—QErp 



X, 



jSSM-non-QE/jj^ 



f SSM ^x) 



(32) 



in Eq. (1281) . This modeling, which we will call SSM-inel in the following, is thus based 
on the assumption, suggested by the RFG, that there exists only a single universal scaling 
function and that the latter can be identified with the phenomenological fit obtained from the 
QE longitudinal data. Henceforth, for simplicity we denote the phenomenological (super-) 
universal scaling function to be used both for impulsive QE and inelastic contributions by 

jSSM 

Important ingredients of the model are, of course, the single-nucleon structure functions. 
In our past work [6j, which was focused on the highly- inelastic scattering region, we used the 
Bodek et al. parametrizations of the proton and neutron structure functions which were 



available at the time. However, in recent years new studies, both theoretical 
experimental [46, 



47 



48 



44 



451 ] and 



491 ] , of the nucleon structure functions in the resonance region have 



been performed, indicating the need for more sophisticated parametrizations. As we are now 
studying this region, we have updated our calculations using more modern expressions for 



"f'2. We thus use parametrizations recently obtained by Bosted and Christy j50|, both for 



w 



the proton [46( and neutron 



471 ] structure functions. We note that all details regarding the 



nucleon resonances, such as finite widths, are automatically included in the parametrization 
and that we do not consider any possible medium-modification of single-hadron properties. 

In order to understand better the role played by higher resonances, we also consider 
a variant of the full SSM-inel model. For this approach, denoted SSM-A, we consider 
contributions coming only from N — > A excitations, which we describe in terms of form 



factors, and we take the finite width of the A explicitly into account. Following [32) we start 
with 



Rf 



L,T 



M2 1 



r(/j,x)/2m 



N 



-R£ T (K,\,ii x )dLi 



x 



/<Ar • r(// v ) 2 , im L v ' ' " w ' 1 (33) 

where R£ t (k, X, fix) are the RFG response functions of Eq. (lT4"j) calculated using a generic 
nucleon excitation invariant mass fix, and tp x is obtained from Eq. (ITS]) for fi& — > fix- Once 
again we then generalize the RFG model by substituting for the RFG scaling function in 
Eq. (fl4l) the universal one, f A (ipx) f SSM {'4 ) x)- The integration limits in Eq. (133!) are 
those of Eqs. (13T1) . and the fix dependence of the A width Y is given by 
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FIG. 3: Cross sections in nb/sr/MeV (left-hand panels; SSM-QE results also included) and non-QE 
superscaling functions (right-hand panels) for 12 C, calculated within the SSM-inel (dashed line) 
and SSM-A (dotted) approaches, and compared with the phenomenological SuSA results (full 



The kinematics selected here are summarized in Table I; the data are taken from from 



34, 



35 



mc 



15 



/fA / p; 



with T = 120 MeV, 



(// 



/4 



(34) 



(35) 



and where p^ es is obtained from Eq. (1551) with /ix = A*a- We then compute inclusive 
cross sections using the response function in Eq. (133]) and, in order to obtain superscaling 
functions within this model, namely f SSM ~ A (ipA), as usual we divide the cross sections by 
S A /kp, where S A is the factor given in Eq. (126]) . These superscaling functions may then be 
compared with the phenomenological SuSA one, fsuSA~ mel {' l P^)^ discussed above [Eq. f[2"5"]) ]. 



B. Results 



In this section we illustrate the results for the superscaling function obtained using the 
SSM-inel and SSM-A models, together with the phenomenological SuSA fit. Before studying 
the behavior of the non-QE scaling function over the whole range of kinematics considered 
in Fig. [21 we will present a few selected examples of cross sections and scaling functions. 
The use of cross sections allows a direct comparison with "real" and more familiar data, and 
the selection of fixed kinematics can illustrate better the characteristics, and the limits, of 
the models. This comparison is shown in Fig. [31 where the left-hand panels show results 
for cross sections and the right-hand panels show results for non-QE scaling functions. For 
illustration, we choose to consider kinematics covering a limited range of energy and mo- 
mentum transfer, large enough so that the A excitation is clearly present and small enough 
so that higher inelastic contributions do not overlap completely with the QE and A peaks. 
Specifically, we select a lower limit case (panels al, a2) corresponding to incident energy 
e = 620 MeV and scattering angle 9 = 36°, and an upper limit case (panels dl, d2) with 
e = 3595 MeV and 9 = 16°. In order to explore the angle dependence of the cross sections 
and scaling functions, in the middle panels we show results for intermediate kinematics with 
two choices of scattering angle, e = 680 MeV, 9 = 60° (panels bl,b2) and e = 1299 MeV, 
9 = 37.5° (panels cl,c2). The kinematics are summarized in Table I which contains as well 
the momentum transfers at the QE and A peaks, qQE and q&, respectively. 

In the figure we compare the results obtained using our SSM-inel and SSM-A models with 
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Case 


e [MeV] 


9 [deg] 


Qqe [MeV/c] 


q A [MeV/c] 


a 


620 


36 


366 


460 


b 


680 


60 


606 


600 


c 


1299 


37.5 


791 


850 


d 


3595 


16.02 


1056 


1189 



TABLE I: C(e,e') kinematics considered 



those corresponding to the SuSA fit introduced at the end of Sec. Ill B[ The SSM-A model 
is certainly an overly simple one and, as can be seen from Fig. [31 the corresponding curves 
show the largest discrepancies with the data. However, the SSM-A results are qualitatively 
interesting because, when compared with the SSM-inel results, they allow us to some extent 
to disentangle the effects related to contributions arising from higher resonances, which 
cannot be eliminated from the data. 

The cross sections plotted in the left-hand column of the figure include the QE contribu- 
tion calculated within the SSM-QE modeling outlined in Sec. Ill which is the same for all 
models. Differences between the various curves in the QE region are therefore due to differ- 
ences in the non-QE part of the cross sections obtained using the various models. By looking 
at the cross sections we can clearly see that our inelastic model always underestimates the 
data in both QE and, especially, inelastic regions. More specifically, for small incident energy 
(upper panels) the QE peak is well reproduced. At the A peak both SSM models clearly 
underestimate the data, while SuSA obviously reproduces the peak reasonably, since it was 
fit to the data. All models are unable to reproduce the cross section completely in the region 
between the QE and A peaks. 

Similar results hold for the A-peak region at larger scattering angles (panel bl). We 
notice that in this case the SSM and SuSA modeling underestimates the data even at 



the QE peak. This is related to the fact that at 



arge scattering angles the transverse 



contribution is dominant. Previous scaling studies [3|, |4( in fact showed that the transverse 
QE superscaling function extracted from the data differs from the longitudinal one and 
exhibits stronger scaling violations. We attribute these differences to contributions beyond 
the impulse approximation, such as 2p-2h MEC and correlations, which are not included in 
the models discussed in this section (see, however, Sec. IIVj) . Moreover, at larger angles the 
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overlap between the QE and A peaks becomes more significant, which explains the difference 
between SuSA and SSM models at the QE peak. 

The same considerations can be extended to the case of higher incident energies (panels 
cl and dl). We observe that in these cases the SSM- A results decrease very rapidly at 
large energy transfer, as does the SuSA curve, because no higher inelastic contributions 
beyond the A are included. The SSM-inel curve has an a;-dependence similar to that of the 
data for large energy transfer, suggesting that the single-nucleon inelastic content has been 
correctly implemented in the model. However, the experimental cross sections are again 
underestimated even in the higher inelastic region. 

With these considerations about cross sections in mind, we can now examine the right- 
hand panels of the figure, which show the non-QE scaling functions. We can summarize our 
findings as follows. As already said, the SSM-inel model always underestimates the data. 
This difference, in both size and shape, is particularly relevant for small incident energies 
and, at all kinematics, for relatively large negative values of the scaling variable At 
very low energy (upper panels of the figure) or for ip' A < — 1 this is expected, because in 
these re gion s effects stemming from correlations and 2p-2h MEC can play an important 



role 



is, 



191 ] and they cannot be reproduced by models which assume impulsive, quasi-free 



scattering on bound nucleons. 

In the region —1 < < the theoretical SSM-inel curves still fall below the data, 
but their shape is similar to that displayed by the experimental scaling function. The 
discrepancies are larger below ipA = ~ 0.5, where 2p-2h MEC may still contribute sizably, 
whereas when approaching the A peak the theoretical curves lie closer to the data. The 
conclusion we draw from these observations is that the basic idea of the phenomenological 
superscaling-based model (SSM-inel) is probably correct and that it can account for most 
of the difference in shape between the experimental QE and non-QE scaling functions, but 
that the model presented here is still too simple and needs some improvements in order to be 
considered quantitatively reliable. In particular, as previously observed, the model assumes 
universality of the longitudinal and transverse QE scaling function, i.e., the so-called scaling 
of the 0th kind, which has been shown to be violated by the QE data. While part of this 
violation can be ascribed to correlation and 2p-2h MEC effects, as discussed in the next 
section, a certain amount of it could be present even at the impulse approximation level, 
and, if so, should be incorporated in the model by using different scaling functions for the 
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T and L responses. This would lead to a renormalization of the calculated cross sections 
and non-QE scaling functions, which may fill some of the discrepancy with the data at the 
A peak. Unfortunately, such an improvement of the model is not straightforward, although 
work is now in progress along this line. 

If we accept that at least some of the difference in normalization between the data and 
the calculated f non -Q E close to the A peak can be accounted for by an improvement in the 
scaling functions used as ingredients in the model, then the SSM-inel results obtained so far 
can provide some useful additional insight on the behavior of the A superscaling function. In 
Fig. H] we plot the function f non -Q E for a relatively large set of kinematics (indicated in the 
key inside the figure) corresponding approximately to values of the momentum transfer in the 
range 500-1500 MeV/c. The top panel shows the experimental non-QE scaling functions, 
the middle panel those obtained within the SSM-inel model and the bottom panel those 
calculated with the SSM-A model. We see that both the data and the SSM-inel scaling 
functions present the same type and degree of scaling violations in the region — 1 < < 0, 
with the curves corresponding to the highest momentum transfer being the lowest ones for 
approximately < —0.5 and then becoming the highest one for larger values of the scaling 
variable. In contrast, this behavior is practically absent in the SSM-A results, suggesting that 
scaling violations in the region — 1 < ip' A < are essentially due to contributions from higher 
resonances. This observation has important consequences for SuSA modeling of neutrino 
cross sections [zj], because it supports the validity of using the universal scaling function 
jSSM - n p rec !i c ting cross sections for kinematical conditions in which only contributions up 
to the excitation of the A resonance are relevant. 

Still looking at Fig. HJ let us mention that both SSM models provide scaling-violations 
at the A peak which seem to be larger than those exhibited by the data. In our study we 
have checked that this is due to kinematical effects, being related to the interplay between 
integration limits and the dependence of the variable ip' x upon the invariant mass fix- The 
inclusion of some degree of scaling-violation in the phenomenological scaling function used 
in the model may solve this problem. 

The differences in the behavior of the theoretical and experimental scaling functions at 
the peak of the A may also be related to the different role of final-state interactions (FSI) for 
QE scattering and A excitation. In fact, previous studies in the QE region have shown that 
the phenomenological QE scaling function is affected by FSI at the right of the QE peak 
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FIG. 4: (Color online) Experimental "QE subtracted" data for f non -Q E for 12 C for a variety of 
kinematics (top panel) and corresponding results of the SSM-inel model (middle) and of the SSM- 
A model (bottom). The kinematics considered are labeled with e (MeV) and 6 (deg). The SuSA fit 
is shown by the solid (red) line. The values of the momentum transfer for the kinematics presented 
here fall approximately in the interval 0.5 < q < 1.5 GeV/c. 
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FIG. 5: Non-QE superscaling function for C calculated within the SSM-inel model versus the 
A scaling variable t/a- The same kinematical cuts as in Fig. [2] are considered, as indicated in the 
different panels. The phenomenological fit of the non-QE superscaling function is also shown for 
comparison (solid line). 

where FSI produce a tail, and partially at the peak, since a larger tail at positive ip results 
in a smaller maximum value of the scaling function. While the tail of the phenomenological 
function contributes very little to the calculated non-QE scaling function at the left of the 
A peak due to the limits of integration (see Eqs. fl28l [3T1 and |33|) ) . its maximum value may 
have some relevance at the A peak. Some details concerning the limits of integration and 
the role of f SSM in determining the non-QE scaling function can be found in the Appendix. 



Finally, before proceeding in the next section to the analysis of the residual after impulsive 
contributions have been removed, and to complete the overview of our results for the function 
jnon-QE ^ we s j low m pjg_ [5] ^he complete set of SSM-inel results for all kinematics where 
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FIG. 6: Full quasielastic scaling function f® E for 12 C as a function of i^'qe- ^he same kinemat- 
ical cuts as in Fig. [2] are considered, as indicated in the different panels. The solid line is the 
phenomenological QE fit used in this work. 

data are available, with the same kinematical cuts used for the results presented in Fig. [2] 
(see also Table I). Also, for comparison, in Fig. [6] we show f QE (Eq. dTDJ) function of 

IV. RESIDUAL NON-IMPULSIVE CONTRIBUTIONS AND SYNTHESIS OF 
THE CROSS SECTION 

The superscaling-based model developed in the previous sections allows one to study the 
behavior of the superscaling function within the context of the impulse approximation and 
therefore to assess the size of any potential non-impulsive contributions. In particular, it 
is interesting to combine the two impulsive contributions denoted SSM-QE (Sec. Ill Ap and 
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FIG. 7: (Color online) Cross sections in nb/sr/MeV versus u. The kinematics corresponding to 
the various labels are the same as in Fig. [3] and Table I. The curves are discussed in the text. 



SSM-inel (Sec. IIII A|) and subtract this from the data to yield a residual: 

exp / j \ SSM-QE / , \ SSM-inel 



da 



da 



da 



da 



\de'dttj \de'dttj \de'dVt) \de'dVt) ' ^ 

The results are shown in Fig. [7] for the same kinematics considered in Fig. [3j Here the 
(black) stars are the complete experimental data, the (red) squares the QE-subtracted cross 
sections [Eq. ( TT3!) ] and the (blue) circles the residual cross sections [Eq. ( 136]) ]. Note that, 
lacking any means of evaluating what errors are incurred in the subtraction procedures, we 
have not given any uncertainties for the non-QE and residual cross sections shown in the 
figure. 

Focusing on the residual cross sections we see significant contributions left over after the 
SSM-QE and SSM-inel results have been removed. As stated several times in the previous 
sections we expect there to be non-impulsive effects from 2p-2h MEC [bsl . \\\ . Indeed, when 
these are compared with the residuals (shown as solid magenta curves in the figure) one 
sees rough agreement. That is not to say that one now has a fully satisfactory picture of 
inclusive electron scattering in this kinematic region — there are still several open issues. 
In particular, when MEC effects are included (and they are not optional; they must be 
included) gauge invariance requires that corresponding correlation contributions must also 
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occur. In 
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271 . l28l . l29l . |30| this problem was dealt with for the lp-lh sector. However, this 



has not yet been done for the 2p-2h response, although work is in progress [5l| to address 
this issue. Another issue goes back to comments made in the previous sections, namely, even 
the SSM-QE approach has some uncertainties in that scaling of the Oth kind may be broken 
to a small degree, and that the somewhat larger transverse scaling functions found in the 
RMF approach and the lp-lh MEC contributions (which lead to a small reduction of the 
transverse cross section) may not completely compensate one another. In effect one could 
break the Oth-kind scaling by using a slightly different scaling function for the transverse 
contributions and thereby modify the residuals seen in the figure. It is clear, however, that 
a significant amount of the residual can be explained by the 2p-2h MEC contributions. In 
the last section we shall return to this point and comment on the implications this has for 
predicting neutrino reaction cross sections. 

In this section we focus primarily on the cross sections and only at the end of the section 
we will briefly return to discuss the non-QE scaling function in order to assess the validity 
of SuSA-based models. Note that we should not expect the 2p-2h MEC contributions to 
scale using either type of scaling discussed above, i.e., either the QE type or the A type. In 
fact these contributions have their own characteristic scaling behavior and work in progress 
is aimed at exploring this behavior in the residual data. With these comments in mind, let 
us work in the opposite direction and, rather than analyzing the cross section, attempt to 
synthesize it using the three types of contributions. 

In Fig. [8] we show the net result of adding together the SSM-QE, the SSM-inel and the 
2p-2h MEC contributions for comparison with the data. The results are quite encouraging: 
the basic qualitative structure of the data is also present in the net result of the superscaling 
analysis, although clearly there is more to be done before one can claim to have a fully 
quantitative description of inclusive electron scattering in this region of kinematics. In 
particular, the net result of adding the three contributions falls short of the data in the QE 
peak region, and this might be fixed by slightly breaking the Oth-kind scaling (as discussed 
above) or by exploiting the flexibility that is inevitably present in the modeling of the 2p-2h 
MEC contributions (for instance, by using a different shift energy than the one that was 
chosen for the results presented here). It should be stressed that this rather good level 
of agreement between theory and experiment has been obtained by adding together three 
separate contributions, each with its own distinctive kinematic dependence, and thus any 
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FIG. 8: Cross sections in nb/sr/MeV versus uj. The kinematics corresponding to the various labels 
are the same as in Fig. [3] and Table I. The curves are the sum of the SSM-QE, SSM-inel and 2p-2h 
MEC contributions. 

attempt to represent experimental data using only a subset of the contributions is bound to 
fail for some choice of kinematics. 

In order to make contact with the Superscaling Analysis of 0], we conclude this section 
by taking the non-QE superscaling functions obtained by using the sum of SSM-inel and 
2p-2h MEC cross sections and inserting them in Eq. ( |25l) . These are shown in Fig. [9J for the 
kinematics considered in the previous figures (see Tabled]). The point of doing this, despite 
the concluding statements made in the preceding paragraph, is to provide a comparison with 
the phenomenological SuSA results discussed in Sec. Ill Bl We observe that the inclusion 
of 2p-2h MEC contributions brings the calculated non-QE scaling function closer to the 
phenomenological fit, supporting the validity of the SuSA-based model for lepton-nucleus 
cross sections at kinematics dominated by A excitation. However, the strength shown by 
the residual data close to the QE peak, not accounted for by the theoretical MEC curves 
considered in this section (as discussed above), affects the non-QE scaling functions at ip' A 
values below approximately -1.5. This can be seen by examining the right-hand column 
of Fig. [3j for instance. These effects should be carefully considered in the future when 
constructing quantitatively reliable models. 
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FIG. 9: Non-QE superscaling functions calculated by using the sum of SSM-inel and 2p-2h MEC 
contributions to the cross sections in Eq. ()25|) . Labels a, b, c, d correspond to the kinematics used 
in Fig. [3] and listed in Tabled The solid line is the SuSA fit. 

V. SUMMARY AND CONCLUSIONS 



In this work we have explored superscaling in electron-nucleus scattering. We have started 
by reviewing the procedures for analyzing scaling in the quasielastic region. A universal 
longitudinal QE scaling function emerges, both based on phenomenology and on modeling. 
Upon assuming that Oth-kind scaling is satisfied (universality of transverse and longitudinal 
scaling functions) we arrive at our model, denoted SSM-QE, for these contributions. 

Next we have focused on the region lying to the right of the QE peak, first introducing 
the definition of the experimental scaling function in this region, f non -Q E . This entails 
subtracting the SSM-QE scaling predictions from the data. We have studied the scaling 
behavior of f non -Q E by analyzing all available high-quality data for 12 C. We have found 
reasonable scaling below the A peak, with scaling violations that can be mainly explained 
in terms of contributions coming from higher resonances. Following the SuSA approach 
presented in jj], we have obtained a phenomenological fit, fs < Hs~A iE i which differs from the 
phenomenological function f SSM ~Q E obtained in previous studies of QE scattering. 

In order to understand this difference, and to explore in detail the breaking of scaling 
shown by f non -Q E ( we have developed an extension for inelastic electron-nucleus scattering 
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within the impulse approximation denoted f SSM ~ inel which is based on previous studies 
of the same type Q|. The model begins with the formulation of the response functions 
in the RFG model, and extends the latter by incorporating in them the universal scaling 
function f SSAI obtained from fits of QE scattering data, making the approach in a sense 
super-universal. The entire inelastic response on the nucleon is incorporated using a recent 
representation of the nucleon's structure functions for kinematics going from pi on-p ro duct ion 



threshold to where DIS takes over 



47]. 



50J, both for the proton [46( and neutron 
The comparison of this impulsive model with the experimental data, for both cross sec- 
tions and non-QE scaling functions, is good, but not entirely satisfactory at first glance, 
since the results always fall below the data. However, the acceptable agreement of the shape 
of the calculated scaling function with the data, suggests that the differences between the 
experimental scaling functions obtained in the QE and A regions could be mainly explained 
in terms of the kinematical effects discussed in the Appendix. Additionally, the results of the 
SSM-inel model allow us to conclude that the scaling violations observed for —1 < ip' A < 
can be mostly explained by the presence of contributions from higher resonances. In par- 
ticular, by comparing results from the full SSM-inel model, in which the entire inelastic 
responses of the nucleons are included, with a variant of this approach (denoted SSM-A), in 
which only the A is included, it has been possible to gain some insight into the roles played 
by excitations lying above the A. 

Having explored the superscaling properties of the SSM-QE/SSM-inel model, we have 
used this model to subtract from the experimental data both the QE contributions and 
those inelastic contributions that can be described within the impulse approximation, thus 
isolating non-impulsive contributions. When this residual is compared with the known 
non-impulsive contributions, namely, those arising from 2p-2h MEC, one sees improved 
agreement between modeling and data. Indeed, it appears that the 2p-2h MEC contributions 
are essential if one is to have a quantitative picture of inclusive electron scattering at the 
kinematics considered in this work. 

Finally, a few words are in order concerning the implications the present study has for 
predictions of neutrino reaction cross sections. Clearly all of the ingredients discussed here 
(QE, inelastic and MEC contributions) also enter in studying the latter, insofar as the 
vector current is concerned. The axial-vector current required for neutrino reactions is 
another matter, however. Unlike the polar-vector current where the leading-order MEC 
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effects enter as transverse effects, but not as longitudinal effects, the axial-vector current 
is the opposite (due to the extra 75 in the basic current, which switches the contributions 
of the upper and lower components in the required matrix elements). Accordingly, for the 
axial-vector currents there are no leading-order transverse effects from MEC, while there 
are for the axial longitudinal/ charge currents. The latter are small for neutrino reactions 
at the kinematics of interest in this work and consequently for neutrino reactions the MEC 
effects enter asymmetrically — essentially, only via the polar-vector currents, but not the 
axial-vector currents. We have seen that the (vector) MEC effects are significant and thus 
any model that does not have them runs the risk of incurring errors of typically 10-20% in 
predicting neutrino cross sections. Using overly simple models such as the RFG is adequate 
for crude estimates of the neutrino cross sections, although almost certainly the SSM analyses 
presented here are considerably better as they capture much of the correct kinematical 
dependences of the polar-vector parts of the electroweak nuclear response. However, as we 
have seen in the present work these impulsive superscaling models do not entirely capture 
all of the necessary content in the currents since they are missing the non-impulsive MEC 
contributions. To the degree that the latter are important one has a complicated problem 
containing at least three parts, the SSM-QE and SSM-inel impulsive contributions together 
with the 2p-2h MEC effects, each with its own distinctive kinematic dependences. 
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APPENDIX: THE A SUPERSCALING FUNCTION AT THE A PEAK 



In Sec. IIII Bl we observed that at the A peak the scaling violations shown by the super- 
scaling function obtained within the SSM-inel model seem to be larger than those present 
in the data. Comparable scale-breaking effects are also obtained within the SSM-A model, 
which includes only the excitation of the A resonance, and therefore they cannot be ex- 
plained in terms of an incorrect treatment of higher resonances. Here we show that the 
origin of these scaling violations in our models is related to kinematical effects and to the 
shape and value of the phenomenological QE function used as input at its peak. In order 
to do so, we work within the SSM-A model, whose simplicity allows us to explore in detail 
the effects of the various terms entering the formulae for the response functions and of the 
corresponding integration limits. 

Let us return to the lower panel of Fig. HJ where we observe excellent scaling for ip' A < 
—0.5, whereas a significant amount of scaling violation remains at larger tp' A and it increases 
approaching the peak. 

In order to understand this behavior let us examine the integral in Eq. together with 
the formulae presented in Sec. Ill Bt We can easily see that, except for minor effects due to a 
residual fix dependence in the coefficients multiplying the form factors in Eqs. ( Hl3f2"0~!) . the 
SSM-A superscaling function / A is essentially given by the integral 

JappxWA) = / -7 ,2 , no Ta 2 J WXjdfXx, (A.1J 

J 1*1 vr - AtA ) + r( / u x ) 2 /4m^ 

whose calculation indeed produces curves that are close to the full f A and have the same 
scaling properties. We can thus use the simpler expression in Eq. ( 1A.1I) to investigate further 
the origin of the scaling behavior of f A and the scaling violations it exhibits at the A peak. 
To do so we choose two fixed values of ipA (for simplicity we also set E s hift = 0), namely 
■0a = (i.e., the Delta peak) and = —0.5 and look at the behavior of the integrand of 
Eq. ( lA.lj) . This integrand is displayed in the upper panels of Figs. [10] and HH for = 
and —0.5, respectively, for different values of the momentum transfer q, as indicated by the 
labels. The asterisk close to origin of the x-axis indicates the integration limit related to 
pion-threshold, while the different dots on the x-axis indicate the upper limit of integration. 
For the largest values of q considered here the latter falls outside the plotted range of fix 
and therefore the corresponding dots do not appear in the figures. 
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FIG. 10: (Color online) Upper panel: integrand of Eq. (|A.1|) . i.e., of fa ppx , as a function of fix, for 
various values of the momentum transfer, as indicated by the labels. The asterisk and points on 
the x-axis indicate the lower and upper limits of integration, respectively. Middle panel: function 
jSSM ag a f unc tion of fix- Lower panel: inelastic scaling variable tpx as a function of fix- All 
curves are calculated at fixed ipA = (here, for simplicity, E s hift has been taken to be zero). 

In the middle panels of the same figures we show f SSM as a function of fix and in the 
lower panels we show how ipx varies with fix for fixed Note that ipx decreases for 
increasing fi X - By looking at Fig. [10] we can see that as q increases, the dependence of tpx 
on fix becomes weaker and thus ipx stays closer to the fixed value of ipA, which in this case is 
0. This means that for larger q the integrand receives contributions mostly from f SSM (ipx) 
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FIG. 11: (Color online) As for Fig. [TU1 but now at fixed V'A = —0.5. 

close to its peak, while for smaller q a more extended range of ipx values contributes. For 
this reason the integrand turns out to be larger for the highest values of q and this gives rise 
to the behavior observed in Fig. HI 

The same type of behavior of ipxi^x) is observed also for ?/>a = —0.5 (lower panel of 
Fig. [TT|) . implying that for high q the variable ipx stays close to the negative value -0.5, thus 
remaining in the region where f s (ipx) is small. In this case, however, this is compensated 
by a larger (with respect to what occurs for smaller values of q) integration interval and thus 
the values of fa Ppx for different q are very close to each other. Moving to even more negative 
values of ipA one could see that the larger integration interval occurring for high q is no 
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longer able to compensate for the smaller values of f SSM (ipx) involved in the integration, 
so that the larger q becomes, the smaller is the value of f£^ px obtained. 

Similar considerations could be applied to the integral in Eq. (1281) entering in the SSM- 
inel model, although in this case a direct comparison between f A and f A ppx cannot be made, 
because of the different single-nucleon ingredients entering the cross section and the dividing 
factor S A of Eq. (|26|) . which defines f A . In particular in this case stronger scaling violations 
appear even in the negative ^-region, which, as discussed in the paper, are related to the 
presence of higher resonance contributions in the inelastic single-nucleon structure functions. 
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